AD-A215  284 


U  NCL  ASSl  FI  ED 


-■  i  ■ 


sic-*  *v  :**ss 


on  of  r  -  s  5  a  :e 


ffr?7 


REPORT  DOCUMENTATION  PAGE 


[TT^omr  security  CLASSi^CATiOs 

"jNCLASS  IF  IED 


*.  slCw«' ’  y  CLASSIFICATION  auThQRiTy 


to  D£C^ASSiF  CATiON/OOWNGAAO»NG  SChEOu^E 


50  aE3TR.CT  vE  MAflUiNOS 


)  O  S  r  a  I  0uT  1  ON- A  v  A>  ^  AQ  u.'  T  Y  OF  «E^O«^ 

APPROVED  FOR  PUBLIC  RELEAS: 
DISTRIBUTION  IS  UNLIMITED 


*£AFC*NiiNG  ORGANISATION  REPORT  NuMSER'Si 


S  MGS'  TOP'  NC  O  BGANiZa  T  ion  »6£0»I  numSERiSI 

APOiiXTR.  8  9'  i  5  K  0 


i*.  NAME  OF  1*6  R»QflgiNG  a  RGAN  .  CA  T  1  ON 

Princeton  Universitv 


pc  OF  C  c  t  3  >  v»8C  v 
■  If  appi.c  30i/  / 


'*  NAME  of  MCN'TOFiNG  ORGANISATION 

AFOSR/NA 


i_  AOORESS  (C»ly.  St»tt  m*4  /.IP  Cn* Ui 

Princeton,  N.J.  08544 


70  AOQRESS  iC  »fv  ^(Qft  ana  /.IP  Coact 

BUILDING  410 

BOLLING  AF3 ,  DC  20332-6448 


NFW{  OF  F  UNOING-SFONSO«' NC 

organisation  AFOSR/NA 


Bb  OF  F  .C£  SvMflO  L 
If  appiugOir  . 


9  PROCUREMENT  INSTRUMENT  I  OCN  T  l  F  i  CA  T  ION  NUMBER 

AFOSR  Grant  88-0120 


Be  AOORCSS  <Cify.  5<«(«  «nd  HP  Coa*> 

BUILDING  410  - 

BOLLING  AFB ,  DC  20332-6448 


'0  50u»C(  OF  Funding  NOS 


PROGRAM 
6  lEw£nt  NO 

61102F 


M  Ti  7lE  'tneiu At  54Ci.ruy  C.OMihcat ton  ■ 

rj)Phvsical  Models  for  SuDersonic  Turbuler) 

■Baunvssry-K.i-sr-srro-ccnre- 

Alexander  J.  Smits 


. . . . 

FROjECT 

TASK 

*0*K  UNIT 

NO 

NO 

NO 

2307 

A2 

13a  Tyre  OF  REPORT 

’30  T  i  mE  COvEREC 

»4  OA-E  OF  R6  °0  R  T  ,Vf  o  .  0*yf 

13  fag£  COUNT 

FINAL 

paom  2/15/88  -c  6/14/89 

1989  October  17 

30 

It  SUFPUEMtNTAR  *  NOTATION 


IT  CCSATi  COCES 

<’UO  1  G»OU»  1  Su8  OR 

!  i 

If.  ABSTRACT  fCOrtimu*  O*  if  ntc«  inrv  Jrtd  J»  Oi  u/f  *»  ■»u>"9r'' 


see  reverse  side 


DTIC 


ELECTE 
DECO  7 1983 


20.  0‘STRlBuTlON/AvAiWABlLlTY  OF  AlSTRAC” 

21  ABSTRACT  SECURITY  CLASSIFICATION  1 

UFtCLASSlFieO'UNUlMITEO  SX^WE  AS  »*t  0  0  r  C  .  SE  °  S  Z 

UNCLASSIFIED 

i  _  _ 

ia«.  NAME  OF  RtSFONSlitt  INOIVlOUAL 

2  2b  rEc6FH.ON6  NUMBER 

23c  OFFICE  Symbol 

1  Inctudf  A'TFCod#; 

JAMES  M  MCMICHAEL 

202-767-4935 

AFOSR/NA 

00  FORM  1473,  83  APR  soition  o«  i  .an  r:  .$  obsolete  U  N  C  L  A  S  S I  FI  Rn 


tiruaiTy  rt  *44,*inn^A.  *  -  *  '*>*/*‘ 


APOtR  TK.  8  9“  1  580 


1  0 

This  report  is  the  final  report  for  AFOSR  Grant  88-0120,  monitored  by 
Dr.  J.  McMichael .  An  experimental  program  was  carried  out  to  study  the 
detailed  structure  of  supersonic  turbulent  boundary  layers.  The  experiments 
were  designed  to  elucidate  physical  models  and  mechanisms  that  are  particular 
to  compressible  turbulence,  such  as  the  effects  of  compressibility  on  the 
nature  of  the  large-scale  motions,  the  scaling  laws  for  high  Reynolds  number 
supersonic  turbulent  flows,  direct  compressibility  effects  that  cause  the 
exchange  of  turbulence  energy  among  the  vorticity,  entropy  and  sound  modes, 
and  the  transport  of  heat  and  momentum  by  compressible  turbulent  motions. 

In  Section  2.  a  detailed  description  the  completed  work  is  given, 
including  a  description  of  the  new  optical  experimental  tools  developed  in 
conjunction  with  Prof.  Miles.  Appendix  A  contains  a  list  of  the  publications 
acknowledging  this  grant,  and  Appendix  B  gives  a  list  of  the  personnel 
involved  with  the  work,  including  the  degrees  awarded  during  the  period  of 
the  grant . 
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1 .  INTRODUCTION 

This  report  is  the  final  report  for  AFOSR  Grant  88-0120,  monitored  by 
Dr.  J.  McMichael.  An  experimental  program  was  carried  out  to  study  the 
detailed  structure  of  supersonic  turbulent  boundary  layers.  The  experiments 
were  designed  to  elucidate  physical  models  and  mechanisms  that  are  particular 
to  compressible  turbulence,  such  as  the  effects  of  compressibility  on  the 
nature  of  the  large-scale  motions,  the  scaling  laws  for  high  Reynolds  number 
supersonic  turbulent  flows,  direct  compressibility  effects  that  cause  th; 
exchange  of  turbulence  energy  among  the  vorticity,  entropy  and  sound  modes, 
and  the  transport  of  heat  and  momentum  by  compressible  turbulent  motions. 

In  Section  2,  a  detailed  description  of  the  completed  work  is  given, 
including  a  description  of  the  new  optical  experimental  tools  developed  in 
conjunction  with  Prof.  Miles.  Appendix  A  contains  a  list  of  the  publications 
acknowledging  this  grant,  and  Appendix  B  gives  a  list  of  the  personnel 
involved  with  the  work,  including  the  degrees  awarded  during  the  period  of 
the  grant. 

2 .  DETAILED  DESCRIPTION  OF  COMPLETED  WORK 

Here,  we  review  our  work  on  the  study  of  the  structure  of  supersonic 
turbulent  boundary  layers  in  zero  and  adverse  pressure  gradients,  and  the 
flow  of  a  turbulent  boundary  layer  over  curved  surfaces .  A  description  of 
the  new  optical  techniques  developed  by  Prof.  Miles  is  also  given. 

2 . 1  The  structure  of  supersonic  turbulent  boundary  layers 

A  considerable  amount  of  new  insight  into  the  turbulence  structure  of 
supersonic  turbulent  boundary  layers  has  been  obtained,  through  the  detailed 
study  of  representative  boundary  layer  flows.  The  data  include  measurements 
of  <u’>,  <(pu) ’>,  <v'>  and  -^oujv' ,  and  they  have  been  reported  in  great 
detail  in  the  recent  AGARDograph  by  Fernholz,  Smits,  Dussauge  and  Finley 
(1988).  In  addition,  measurements  of  higher  order  moments  have  been  made, 
and  spectral  data,  probability  density  distributions,  and  correlations  have 
also  been  obtained  in  all  cases.  In  this  respect,  these  data  sets  are  the 
most  complete  data  sets  for  zero  and  adverse  pressure  gradient  supersonic 
boundary  ^yers  currently  available. 

During  the  grant  povinH  of  o  r :i.:u  procure  orafH°r;t|  ann 
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adverse  pressure  gradient  flow  on  a  flat  plate  were  completed  (Fernando 
(1988),  Spina  (1988)),  as  well  as  a  study  of  the  behavior  of  a  supersonic 
turbulent  boundary  layer  on  curved  surfaces  (Donovan  (1989)).  This  work  was 
used  to  assess  the  significance  of  variable  density  on  the  structure  of  high 
speed  turbulent  boundary  layers.  For  the  purpose  of  comparison  with  the 
subsonic  case  we  used  the  results  from  Alving  (1988)  because  of  the  wide 
variety  of  data  obtained  in  that  flow.  The  flow  conditions  for  these  cases, 
and  some  others,  are  summarized  in  Table  I.  Because  of  the  experimental 
conditions  of  these  tests,  and  the  lack  of  high-quality,  detailed  data  from 
other  sources,  it  is  not  always  possible  to  make  a  clear  distinction  between 
Mach  and  Reynolds  number  effects. 

According  to  Morkovin  (1962),  the  dynamics  of  compressible  turbulent 
shear  layers  follow  the  incompressible  pattern  closely,  as  long  as  the 
f luctuauating  Mach  number  remains  small.  In  supersonic  Flows  1  and  2,  the 
freestream  Mach  number  is  in  the  range  where  Morkovin' s  hypothesis  should 
apply  over  most  of  the  boundary  layer  thickness .  The  Mach  number  gradient 
near  the  wall,  however,  is  very  high  (for  example,  in  Flow  1  the  sonic  line 
is  located  at  approximately  y  =  0.0056  =  0.13mm  =  75  t/w/uT),  and  it  may  be 
expected  that  the  fluctuating  Mach  number  can  exceed  unity  in  the  region  of 
maximum  turbulence  production,  and  Morkovin' s  hypothesis  may  break  down  in 
the  near-wall  region. 

The  simplest  comparison  between  the  turbulence  behavior  in  subsonic  and 
supersonic  boundary  layers  is  to  compare  the  distributions  of  <u'>.  When 
normalized  by  uT  (=  J rw/pw) ,  the  distributions  appear  to  show  a  strong  Mach 
number  effect  (see,  for  example  Schlichting  1968,  p.  659).  However,  if  the 
results  for  Mach  numbers  less  than  5  are  normalized  by  a  velocity  scale 
derived  using  the  wall  stress  and  the  local  density  (=  Jrw/p),  as  suggested 
by  Morkovin,  the  Mach  number  dependence  is  no  longer  evident.  At  a  Mach  number 

of  6.7,  however,  Owen  et  al .  found  that  this  transformation  did  not  seem  to 
collapse  the  data,  indicating  that  hypersonic  flow  may  display  strong 
compressibility  effects.  It  is  also  known  that  the  flatness  profile  in 
supersonic  flows  is  more  constant  than  in  the  corresponding  subsonic  case  (see 
figure  2).  One  uefinition  of  intermittency  is  I  =  (3/Flatness),  where 
Flatness  =  u’v(u'z)  ,  and  the  results  imply  that  the  intermittency  profile  is 
fuller  in  supersonic  flows.  In  contrast,  the  Rayleigh  scattering  density 


Subsonic 

flow  Supersonic 

flow  1  Supersonic 

flow  2  Supersonic 

flow  3 

0.1 

2.9 

2.97 

6.7 

Uref 

31  m/s 

565  m/s 

594  m/s 

1110  m/s 

cf 

.00283 

.00114 

.00166 

.00080 

R* 

50uU 

80,000 

15000 

8500 

6 

19.3  mm 

28  mm 

12  mm 

33  mm 

Table  1 .  Experimental  conditions  for  the  supersonic  boundary  layers 
investigated  by  Spina  (1988)  and  Fernando  (1988)  (Flow  1),  Robinson  (1986) 
(Flow  2),  and  Owen  et  al .  (1975)  (Flow  3),  and  the  subsonic  boundary  layer 
investigated  by  Alving  (1988) . 

field  flow  visualization  by  Smith  at  Princeton  (private  communication) 
reveals  that  the  boundary  layer  has  a  very  similar  appearance  to  that  seen  in 
subsonic  flows  (see  figure  3) ,  and  it  may  be  possible  that  the  flatness 
profiles  of  (p u) '  do  not  give  an  accurate  picture  of  the  intermittency . 

Recent  measurements  of  the  intermittency  function  from  the  (p u) '  signal  by 
Selig  (1988)  seems  to  support  this  idea.  What  is  even  more  interesting  is 
that  the  shear  correlation  coefficient  Ruv  (=  -u'v’  /<u'Xv'>)  is  different 
(figure  4).  The  subsonic  data  reveal  a  higher  correlation  across  the 
boundary  layer  with  a  nearly  constant  value  of  0.45  for  0  <  y/6  <  0.8,  while 
the  supersonic  correlation  decreases  steadily  as  y/6  increases.  Differences 
were  also  observed  in  the  distribution  of  the  structure  parameter  a^  -  - 

-1  O 

u'v'/q  (not  shown).  In  these  same  studies,  however,  Fernando  and  Alving 
found  similar  anisotropy  ratios,  suggesting  that  the  difference  in  character 
of  the  Ruv  distribution  is  caused  by  a  change  in  the  shear  stress,  that  is, 
the  organized  motions,  not  by  a  change  in  u'  or  v'  alone.  These  changes  are 
clearly  evident  in  the  joint  probability  density  distributions  of  u'  (or 

L  n  ■ 

(pu) ’ )  and  v' ,  shown  in  figure  5.  The  contributions  to  -u'v'  are  organized 
differently;  as  a  basis  for  comparison,  it  may  be  seen  that  the  major  axes  of 


these  approximately  elliptical  distributions  are  aligned  more  closely  with 
the  horizontal  axis  when  the  flow  is  supersonic.  A  wide  variety  of  two-point 
space-time  correlation  data  is  available,  including  mass-flux  (pu) '  (or 
velocity  u')  correlations  in  the  streamwise,  normal,  and  spanwise  directions 
(x.y,  and  z,  respectively).  To  begin  the  discussion  of  the  results, 
consider  the  time  records  of  (p u) '  obtained  in  the  zero  pressure  gradient 
boundary  layer  from  three  hot  wires  (figure  6).  The  signals  exhibit  a  very 
similar  character,  indicating  the  passage  of  organized  motions  of  a  scale 
larger  than  the  separation  distance  between  the  top  and  bottom  wires. 

The  space -time  correlations  for  one  supersonic  and  one  subsonic  flow  are 
shown  in  figure  7.  For  both  cases,  the  peak  values  of  the  correlations  are 
quite  high,  reaching  a  maximum  of  0.65  near  the  middle  of  the  boundary  layer. 
The  correlation  frictions  for  the  supersonic  boundary  layer  are  considerably 
narrower  than  for  subsonic  layer.  Furthermore,  the  dimensionless  delay  time 
corresponding  to  the  peak  of  the  space-time  correlation,  rmax ,  decreases  from 
0.4  (  -  20  +  0.5pts)  at  the  floor  to  nearly  zero  at  the  edge  of  the  boundary 

layer . 

The  high  peak  level  of  the  correlation  and  the  non-zero  value  of  the 
time  delay  imply  that  both  wires  are  detecting  the  same  "disturbance",  and 
that  one  wire  is  detecting  it  before  the  other.  Since  the  time  shift  was 
applied  to  the  upper  wire,  the  peak  at  negative  time  delay  means  that  the 
upper  wire  detects  the  disturbance  first,  that  is,  the  disturbance  leans 
downstream.  Accordingly,  an  angle  6  can  be  defined  for  this  "front"  by  using 
the  value  of  rmax  along  with  the  wire  separation  distance,  £,  and  the  local 
convection  velocity.  That  is. 


=  tan 


-1 


[ 


± 


^crmax 


] 


The  angle  6  may  be  called  an  "average  structure  angle,"  in  that  it  is  asso¬ 
ciated  with  an  average  large-scale  motion.  Figures  8  and  9  show  that  the 
structure  angle  depends  on  the  distance  between  the  two  measurement  points, 
and  that  the  distribution  for  supersonic  flow  is  different  from  that  in 
subsonic  flow.  In  the  supersonic  case  for  small  values  of  £  (£/<5  -  0.09 
, say),  the  structure  angle  is  approximately  constant  at  a  value  between  45° 
and  50°  for  0.2  <  y/5  <  0.8.  For  f/8  >0.2,  the  structure  angle  becomes 


insensitive  to  variation  in  the  separation  distance,  and  it  varies  from  about 
40°  at  y/8  =  0.2  to  about  60°  at  y /&  =0.8.  In  the  subsonic  case,  however, 
small  values  of  £  give  larger  values  of  9,  and  large  values  of  £  gives 
smaller  values  of  8  than  those  observed  in  the  supersonic  case.  For  example, 
Alving  found  that  tor  £/5  =  0.1,  the  structure  angle  was  approximately  at  a 
value  of  60°  for  0.2  <  y/8  <  0.8,  whereas  for  f/6  >0.2,  J  varied  from  about 
20-30°  to  about  50°  over  the  same  interval  in  y.  Note  that  the  uncertainty  in 
the  convection  velocity  will  not  affect  these  results  significantly:  when  8 
=  95°,  0%  error  in  Uc  leads  to  an  error  in  6  of  only  2-3°.  These  results 
help  to  explain  the  differences  seen  in  the  correlation  functions  (figure  7). 
Since  the  average  large-scale  structure  is  more  upright  in  a  supersonic 
boundary  layer,  the  space-time  correlation  is  apt  to  fall  off  at  shorter  time 
delays  since  the  extent  of  the  structures  in  the  mean  flow  direction  will  be 
smaller . 

Using  the  same  upstream  flow  conditions  used  by  Spina  and  Smits  for  the 
zero  pressure  gradient  layer  discussed  above,  Donovan  and  Smits  (1987) 
investigated  the  mean  structure  angle  distribution  following  a  short  region 
of  concave  surface  curvature  using  two  different  flow  models:  one  which 
turns  the  flow  through  t>°  with  a  radius  of  curvature  of  1270  mm  (6/R  =  .022), 
and  the  other  turns  the  flow  16°  with  S/R  =  .08.  Fernando  and  Smits  (1987) 
made  similar  measurements  on  a  flat  plate  following  a  short  region  of  adverse 
pressure  gradient.  In  that  case,  the  pressure  gradient  was  generated  by  a 
contoured  plate,  designed  so  that  the  pressure  distribution  matched  that  of 
the  8°  model  using  by  Donovan  and  Smits.  The  general  shape  of  the 
distribution  remained  the  same  as  in  the  zero  pressure  gradient  case. 

However,  there  appeared  to  be  a  small  increase  in  the  structure  angle  after 
each  of  the  three  perturbations.  Furthermore,  the  structure  angle  after  the 
stronger  curvature  was  slightly  higher  than  after  the  weaker  curvature. 
Donovan  and  Smits  suggested  that  the  perturbation  rate  was  too  rapid  in  the 
stronger  curvature  model  to  allow  read’  -stment  of  the  large-scale  motions, 
and  thus  the  angle  of  inclination  is  affected.  It  appears  that  the  same 
preliminary  conclusion  can  be  drawn  for  all  three  of  the  flow  perturbations, 
since  they  all  exhibit  the  same  trend. 

While  the  present  study  traversed  two  "detection  probes"  through  the 
boundary  layc^r  at  a  fixed  separation  distance  (small  compared  to  8),  most 
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other  measurements  of  this  kind  have  used  one  detention  probe  fixed  at  the 
wall  (a  hot  wire,  a  shear  stress  gauge,  or  some  similar  device)  and  another 
probe  which  was  traversed  through  the  boundary  layer,  thereby  varying  the 
separation  distance.  The  fixed  separation  method  used  here  results  in  a 
tvpical  mean  structure  angle  of  85°  in  supersonic  flow  and  about  30°  in 
subsonic  flow.  While  the  variable  separation  method  seems  to  give  a  lower 
characteristic  value;  in  superscni''  flow  Robinson  (1988)  found  30°  and  in 
incompressible  flows  Brown  and  Thomas  (1977)  found  18°,  whereas  Rajagopalan 
and  Antonia  (1979)  found  12.5°,  and  Robinson  (1985)  found  16°.  The  advantage 
of  the  present  method  is  that  the  slope  of  the  structure  is  determined 
locally,  instead  of  being  inferred  from  a  large-scale  measurement. 

A  superposition  of  peak  time  delays  on  the  "mean  structure  shape"  was 
used  to  produce  iso-ccrrelation  contours  for  the  supersonic  case  using  f  = 
0.096  and  0.306  (figure  10).  The  "mean  structure  shape”,  the  solid  curve 
drawn  through  the  center  of  the  contours,  was  constructed  from  an 
extrapolation  of  the  mean  structure  angle  from  one  measurement  location  to 
the  next  /"with  the  first  location  supplying  the  appropriate  inclination  from 
the  origin).  The  peak  of  the  cross -correlation  at  each  mean  wire  position 
was  then  shifted  so  that  it  was  coincident  with  the  curve  delineating  the 
mean  structure  shape.  The  contours  give  a  good  indication  of  the  extent  of 
the  field  around  the  identified  structure.  The  same  mean  structural 
characteristics  are  evident  in  ail  of  the  contours,  as  well  as  those 
determined  by  Robinson.  Similar  contour  plots  were  derived  for  the  subsonic 
boundary  layer  studied  by  Alving  and  they  are  given  in  figure  11.  The 
differences  in  mean  structure  angle  that  exist  between  supersonic  and 
subsonic  flows  (figures  8  and  9)  are  readily  apparent  in  figures  10  and  11, 
as  are  the  much  greater  streamwise  extent  of  the  large-s  ale  structures  in 
the  subsonic  boundary  layer. 

Donovan  f 1989)  has  studied  the  space-time  correlations  in  the  supersonic 
boundary  layer  flow  downstream  of  the  short  region  of  concave  surface 
curvature  described  earlier.  He  observed  a  significant  elongation  of  the 
average  structure  in  the  streamwise  direction,  consistent  with  the 
distortion  that  would  be  observed  if  the  large  scale  motion  can  be 
represented  by  a  horseshoe  vortx  model. 

In  addition  to  the  measurements  made  at  two  points  separated  in  the 
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direction  norm.il  to  the  wall,  measurements  were  taken  at  a  variety  of 
spanvise  spacings  in  both  the  supersonic  and  subsonic  cases.  The  space -time 
correlation  for  a  spanwise  spacing  of  0.096  is  shown  in  figrre  12  for  several 
locations  across  the  supersonic  boundary  layer.  The  corresponding  results 
for  the  subsonic  boundary  laver  are  given  in  figure  13.  The  peak  values  are 
somewhat  lower  than  in  the  supersonic  case,  although  the  correlations  are 
much  broader.  Tn  both  cases,  however,  the  character  and  strength  of  the 
correlations  is  similar  to  those  found  for  the  vertical  separations,  except 
that  the  peak  occurs  at  zero  time  delay  in  this  case.  Since  the  correlation 
between  wires  spaced  0.096  apart  is  similar  for  both  the  spanwise  and 
vertical  alignments,  a  similarity  of  the  structures  is  suggested  in  the  y- 
and  z-directions  lor  small  distances  in  both  cases.  As  determined  from  the 
peak  values  of  the  correlation  functions,  it  appears  that  the  spanwise  scales 
arc  slightly  smaller  than  the  vertical  scales. 

A  comparison  of  the  length  scales  based  entirely  upon  the  peak  values  of 
t In  correlations  is  not  conclusive,  however.  Therefore,  the  transverse 
scales  were  further  explored  with  the  aid  of  iso-correlation  contours  which 
give  the  behavior  of  the  entire  correlation  curve,  not  just  the  peak  value. 
The  space -time  correlations  were  reflected  about  £  =  0  for  each  y-position 
1  this  is  valid  since  the  spanwise  correlations  are  symmetric),  anu  iso¬ 
correlation  contours  were  drawn  from  the  resulting  surfaces  (figures  19  and 
13;.  As  with  the  vertical  correlat ions ,  the  time  delay  was  normalized  by 
outer- laver  variables,  and  by  Taylor's  hypothesis  can  be  interpreted  as  a 
streamwise  distance.  The-  plots  then  give  (pseudo)  x-z  cross-sections  of  the 
boundary  laver  at  three  different  y- locations.  What  is  very  striking  is 
that  the  spanwise  extent  of  the  large-scale  motions  in  subsonic  and 
supersonic  flows  are  almost  identical  (and  in  good  agreement  with  the  results 
of  Kovasznav  et  al .  19/0),  whereas  the  streamwise  scale  differs  by  a  factor 
of  about  two.  In  addition,  the  spanwise  scale  of  the  detected  structures 
increases  away  from  the  wall.  Since  the  overall  size  of  the  structures 
increases  with  y/6  this  behavior  is  not  surprising.  The  spanwise  scale  of 
the  detected  organized  structures  should  therefore  increase  as  we  move 
farther  from  the  wall. 


All  the  measurements  presented  here  indicate  that  despite  broad 
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similarities,  the  turbulence  structure  of  supersonic  and  subsonic  boundary 
layers  display  significant  differences.  Some  of  these  differences,  such  as 
the  change  in  the  flatness  profile,  have  been  observed  in  previous  studies, 
and  are  now  relatively  well  known.  Structure  parameters  have  not  been  widely 
studied  in  supersonic  flows,  however,  and  the  new  measurements  presented  here 
indicate  that  strong  differences  exist.  For  example,  the  length  scales 
derived  from  space-time  correlations  indicate  that  the  spanwise  scales  are 
almost  identical  but  that  the  streamwise  scales  in  the  subsonic  flow  are 
about  half  the  size  of  those  in  supersonic  flow.  The  large-scale  structures 
in  the  subsonic  boundary  layer  also  appear  to  move  slightly  slower,  and  lean 
more  towards  the  wall,  than  those  observed  in  supersonic  flows,  and  their 
shear  stress  content  is  distributed  differently  among  the  four  quadrants. 

All  these  observations  suggest  that  there  may  be  fundamental  differences 
between  the  structure  oi  subsonic  and  supersonic  boundary  lavers .  It  is 
possible  that  the  density  gradients  in  a  supersonic  shear  layer  affect  the 
large  scale  structure,  and  that  there  exists  a  damping  effect  of  Mach  number 
on  the  turbulent  motions  which  may  be  important  even  for  turbulence  away  from 
the  wall.  In  some  sense  it  would  not  be  surprising  to  find  differences 
between  compressible  and  incompressible  boundary  layers  since  the  vorticity 
transport  equation  describes  the  transport  of  vcrticity  per  unit  mass,  father 
than  the  absolute  vorticity,  and  there  is  a  contribution  from  the  non- 
barotropic  term.  Density  gradients  must  therefore  affect  the  vorticity 
dynamics  to  some  extent,  and  the  extent  of  this  influence  will  vary  with  Mach 
nuinbe  r  , 


2.1  New  Optical  Techniques  For  Investigating  Compressible  Turbulent  Flows 
The  image  of  the  shock  wave/boundary  layer  interaction  shown  in  figure  16  was 
made-  possible  by  using  a  high-power,  far  ultraviolet  excimer  laser  in 
conjunction  with  a  high  -  sens itivity ,  far  ultraviolet  camera.  By  focusing  the 
laser  into  a  thin  sheet  of  light  and  passing  it  through  the  wind  tunnel, 
cross-sectional  images  of  the  air  density  can  be  recorded  by  direct  Rayleigh 
scattering.  Due  to  the  tv +  dependence  of  the  Rayleigh  scattering  cross 
section,  operation  in  the  far-UV  produces  a  signal  which  is  approximately  60 
times  that  in  the  visible.  Even  though  the  excitation  is  near  resonant 
transitions  in  oxygen,  the  Rayleigh  scattering  cross  section  of  air  in  this 
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region  is  virtually  independent  of  laser  frequency  and  molecular  temperature 
(Miles  et  al  .  1988a).  As  a  consequence,  the  collected  light  level  is  a 
direct  measure  of  the  air  density.  Illumination  is  with  an  argon- fluoride 
laser  operating  in  the  vicinity  of  .193  microns  with  a  pulse  duration  of  10 
nsec.,  so  the  cross-sectional  image  is  frozen  in  time.  This  diagnostic 
technique  mav  be  used  together  with  vibrational  tagging  of  oxygen  molecules 
to  yield,  simultaneously,  an  instantaneous  velocity  profile.  This  method  for 
generating  velocity  profiles  combines  Raman  Excitation  plus  Laser - Induced 
Electronic  Fluorescence  (RELIEF)  and  some  recent  results  for  experiments  in 
free  shear  layers  were  reported  by  Miles  et  al .  (1988b).  Examples  of  lines 
which  have  been  marked  before  and  after  the  Mach  disk  of  an  underexpanded 
fret,  jet  are  shown  in  figure  17.  The  lines  in  both  cases  were  marked  2  /jsec 
be  tore  the  image  was  taken,  so  the  line  deformation  is  a  quantitative  measure 
of  the  velocity  profile.  Rayleigh  scattering  was  simultaneously  recorded  to 
give  the  density  cross  section.  The  Rayleigh  scattering  technique  has 
already  given  the  first  detailed,  quantitative,  two-dimensional  picture  of 
the  turbulent  structure  of  a  supersonic  boundary  layer  (see  below).  The 
RELIEF  method  is  currently  being  extended  to  boundary  layer  measurements  to 
obtain  simultaneous  time-lines  which  will  complement  the  density  cross 
sections  presented  here. 

A  small  Mach  3  supersonic  wind  tunnel  was  constructed  and  fitted  with  UV 
transmitting  quartz  windows  so  that  the  laser  sheet  could  be  passed  through 
the  flow  field.  The  high-sensitivity  camera  observed  the  scattering  at  90° 
and  images  could  be  recorded  up  to  the  camera  framing  rate  of  30  Hz.  At  the 
point  of  interrogation,  the  Reynolds  number,  based  on  the  boundary  layer 
momentum  thickness,  was  15,000.  The  tunnel  was  operated  as  a  blow  down 
facility  and  exhausted  into  atmospheric  pressure.  An  example  of  a  single 
cross-sectional  image  taken  normal  to  the  boundary  layer  in  the  streamwise 
direction  is  shown  in  figure  18  (the  flow  is  the  same  as  that  shown  in 
figure  3:  the  flow  is  from  right  to  left).  Regions  of  high  density  are 
bright  and  regions  of  low  density  are  dark.  There  is  some  scattering  from 
the  wall  which  causes  a  bright  region  at  the  very  bottom  of  the  image.  The 
scale:,  from  top  to  bottom,  is  7  mm  and  the  resolution  is  on  the  order  of 
several  hundred  microns.  The  boundary  layer  thickness  is  approximately  5  mm. 
Figure-  19  is  a  false  color  rendition  of  this  same  boundary  layer  image  with 
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the  scattering  from  the  bot'"°.r".  trimmed  off.  Cross-sectional  images  parallel 
to  the  boundary  layer  (plan  views)  can  also  be  recorded,  and  figure  20  is  a 
typical  plan  view  image  taken  3  mm  above  the  wall. 

In  a  qualitative  fashion,  these  images  show  that  there  are  well-defined 
individual  structures  located  in  the  boundary  layer.  The  dark  areas  of  the 
image  represent  low  density  regions  characteristic  of  high  temperatures.  The 
bright  areas  are  the  high  densities  found  in  the  free  stream.  Intermittent 
penetration  of  freestream  fluid  deep  into  the  boundary  layer  is  clearly 
visible  in  figure  18.  The  plan  view  (figure  20)  gives  an  even  better 
appreciation  of  the  intermittency  at  a  given  height  in  the  boundary  layer. 
From  many  such  images  a  statistical  picture  of  the  intermittency  can  he 
generated . 

Furthermore,  since  the  brightness  of  the  image  is  a  quantitative  measure 
of  the  density  of  the  gas,  these  images  may  be  examined  to  give  probability 
density  distributions  across  the  flow  and  at  different  height  locations 
within  the  boundary  layer.  Density  correlations  can  also  be  developed  to 
give  a  quantitative  measure  of  the  scale  of  the  turbulent  structures.  Figure 
21  is  the  probability  density  distribution  function  taken  from  figure  20. 

Even  though  it  comes  from  a  single  frame,  this  probability  density 
distribution  of  the  density  fluctuations  is  similar  to  the  probability 
density  distribution  of  the  mass- flux  fluctuations  by  Hayakawa  et  al .  (1984) 
in  a  similar  flow,  giving  indirect  support  for  the  Strong  Reynolds  Analogy. 

In  addition,  we  can  now  contrast  this  statistical  result 
with  an  instantaneous  picture  of  the  density  distribution.  It  is  interesting 
to  note  from  figure  19,  for  example,  that  the  boundary  layer  appears  to 
consist  of  large  regions  of  relatively  uniform  low  density,  separated  from 
the  high  density  region  by  a  rather  thin  buffer  zone. 

In  summary,  we  have  obtained  instantaneous  two-dimensional  density  cross 
sections  in  supersonic  flow  fields  using  a  new  nonintrusive  flow  diagnostic 
method.  These  images  have  already  provided  new  insight  into  the  structure  of 
the  boundary  layers,  and  have  generated  statistical  information  on  the 
density  fluctuations.  These  represent  the  first  direct  measurements  of  the 
density  fluctuations  obtained  in  a  supersonic  turbulent  boundary  layer.  When 
combined  with  the  flow  marking  technique,  we  will  have  the  first  opportunity 
to  directly  observe  the  coupling  between  vorticity  and  density. 
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rig.  1.  Distribution  of  Ha,  u ,  v  /  /  tw/p  in  compressible  turbulent  boundary 
layer  with  freestream  Mach  number  of  2.9,  6  -  26  mm,  and  Rg  -  80,000  (outer 
layer  variable). 
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Fig.  2.  Comparison  of  the  flatness  profiles  for  the  streamwise  fluctuations 
from  supersonic  (Spina  1988)  and  subsonic  (Alving  1988)  turbulent  boundarv 


Fig.  3  Density  cross-section  of  a  Mach  2.5  boundary  layer.  Flow  is  from 
right  to  left. 
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Fig.  ^  Comparison  of  the  shear  correlation 
coefficient  (R^y  -  u' v 1  /<u'Xv'>)  in  eupertonic 
(Femendo  1988)  end  subsonic  (Alvlng  1966) 
turbulent  boundary  layer* 


F1r.6.  Three  simultaneously  measured, 
instantaneous  mass-flux  signals  (Flow  1 
Spina  and  Salts  1987) . 


Fig.  5~-  Probability  density  function  for  (a) 
(eu) 'v'  froa  Fernando  (1988),  and  (b)  u'v'  froa 
Alvlng  (1988). 


fig-  t(«)  Cross-correlation  measurements  la  the 
supersonic  boundary  layer,  {/lag  -  0.09.  Proei 
Spina  (1988). 

Fig.  ^(b)  Cross -correlation  at  several 
positions  In  the  subsonic  boundary  layer, 
f/<99  -  0.10.  Froa  Alvlng  (1988). 


Equi -value  space -tine  correlation 
contours  froc  spanvise  separated  Bass -flux 
signals  In  the  supersonic  boundary  layer 
a'  y/6  -  0.20,  b)  y/6  -  0.51,  c)  y/6  -  0.82 
(Spina  1968 )  . 


Fig./^  Equl  -value  space-tUee  correlation 
contours  froa  apanvlae  separated  velocity  signal 
in  the  subsonic  boundary  layer  studied  by  Alvlng 
a)  y/6  -  0.22,  b)  y/6  -  0.52,  c)  y/6  -  0  80. 
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tht  instantaneous  density  field  for  a  shock  wave 
layer  interaction  produced  by  a  compression  corner. 
ic  cross-section  normal  to  the  wall,  obtained  by 
itr,  using  Rayleigh  scattering.  The  flow  is  from 
left,  the  upstream  freestream  Mach  number  is  2 .5" , 
ir.r  boundary  layer  has  Rr'  =  25,000,  and  the  corner 


Figure  !?  Cross-sectional  inage  of  density  and  line 
velocity  profiles  in  an  underexpanded  air 
jet.  Lines  have  been  tagged  before  and 
after  the  Mach  disk  2  ysec .  before  the 
ir.ace  was  recorded. 


Figure  If  Colc-r-enhau.-ed  image  of  th>.  cross  section  across 

the  hound  :-.rv  lave  r  . 


